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In an earlier paper [1] the advantages of incorporating an 
open-ended or weakly-sealed squeeze-film bearing in a flexible 
support structure simulating an aero-engine assembly were 
examined. In the present paper attention is given to 
empirically modelling the hydrodynamics of the more usual 
tightly-sealed squeeze-film bearing, with a view to assessing 
its damping performance. 


INTRODUCTION 

A squeeze- film is an annulus of oil supplied between the 
outer race of a rolling-element bearing (or the bush of a 
sleeve bearing) and its housing. It is often used as a multi- 
directional damping element for the control of rotor 
vibrations . 

In most aeroengine applications, the purpose of the 
squeeze film is to introduce damping so that the rotor can 
safely negotiate any critical speeds and operate smoothly at a 
higher speed. The amount of damping can be quite critical: If 
it is too little, excessive movement will take place within the 
squeeze-film annulus, if too great the damper will act like a 
rigid link, making its presence superfluous* 

A squeeze- film damper is in effect a journal bearing in 
which the inner member does not rotate. Hence, it is not 
possible to assign to it any linear stiffness coefficients 
which might allow the oil film to take any gravity or other 
static load, since such coefficients are directly dependent on 
journal speed. Any lift emanates instead from non-linear 
effects. For certain operating conditions, such as in the 
regions of critical speeds, where any static load is small 
compared to the dynamic load, circular concentric orbits could 
possibly be assumed, with the result that quasi-linear 
amplitude-dependent coefficients could be used. Linear 
analyses could, thus be carried out to obtain the vibration 
orbits in an interative fashion and transmitted forces could be 
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computed . For cases where the static load is not small a non- 
linear analysis is called for. Any analysis requires a 
specification of the effect of cavitation, which has a profound 
influence on the dynamic load-carrying capacity. Indeed, an 
uncavitated film may be shown to lead to the complete 
elimination of load-carrying capacity. 

In order to investigate closely the performance of the 
squeeze-film damper in this kind of application, an 
experimental rigid rotor [1] was designed and built to simulate 
an aeroengine gas turbine running at speeds up to 6000 rev/min 
in rolling-element bearings, one mounted into a housing with an 
annular gap constituting the squeeze-film damper. This damper 
was provided with a circumferential groove and supplied with 
oil at a prescribed pressure. Each land on either side of the 
groove was of diameter 130 mm and length 9 mm and the radial 
clearance of the squeeze-film was 0,216 mm. The purpose of the 
test rig is to study the orbits of vibration and dynamic 
squeeze-film pressures resulting from various degrees of mass 
unbalance . 


Notation 

c radial clearance of squeeze film 

d seal gap 

e journal eccentricity 

f(z) function of z 

h squeeze-film thickness 

h nondimensional film thickness = h/c = 1 + ecosQ 

l radial width of seal 

L bearing-land length 

m effective rotor mass 

p oil film pressure 

p c ( 9 ) pressure using the long-bearing approximation 

p 0 pressure at maximum film thickness 

p s supply pressure 

Pcav cavitation pressure 

p nondimensional film pressure = p/naj(R/c) 2 

p 0 nondimensional pressure at maximum film thickness 

= p 0 /nu)(R/c) 2 

p s nondimensional supply pressure = p s /r)w(R/c) 2 

Pcav Pcav/ R /c) 2 

P static load per land 

P c dynamic load per land = muoj 2 
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squeeze-film forces in radial, transverse directions 
P/mc<D 2 

P c /mcto 2 = u/c 
bearing radius 
time 

eccentricity of mass unbalance 

axial coordinate 

z/L 

attitude angle measured in direction of unbalance 
force rotation 

nondimensional viscosity = nR(L/c ) 3 /mw 
end leakage flow coefficient 
eccentricity ratio = e/c 
oil viscosity 
circumferential coordinate 
oil film boundaries 
end-leakage factor * (c/d) &/L 

angular velocity of rotor. 

THEORETICAL TREATMENT 


The Reynolds equation, suitably adapted for a squeeze-film 
bearing in which the inner member does not rotate, may be 
written in a nondimensional form as 


Igfh 3 If) + (£) 2 -Wh 3 4E) 

df) dQ L 9 Z 92 


12 ( e 1 cos0 + eot'sine) , 


where 


and 


p = p/nw { §) 2 , 

h = h/c 1+SCOS0 , 

z = z/L , 

( ' ) = d/d ( tut ) . 


( 1 ) 


The variables are given in the notation. 

A common type of squeeze-film damper configuration with a 
central circumferential oil groove and without end seals was 
investigated in ref Cl]. For such a damper, it was assumed 
that the oil flow in the circumferential direction was 
negligible when compared with the flow in the axial direction. 
This is known as the short-bearing approximation. Hence the 
first term on the left-hand side of eqn (1) was dropped. 
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Integrating twice with respect to z, and putting. 


and 

gave 


P = Ps at z = “ 1/2 

p = 0 at z = 1/2 , 


— , n 2 (e'cosQ + sa'sin0) , — 2 1\ — ,1 — \ 

p(6,z) = 6(^) — — -3 (z - -5) + p s ~ z) (2) 


(1 + ECOS 0 ) 


This is the general expression for nondimensional pressure in a 
short squeeze-film with a circumferential oil supply groove and 
without end seals . 


The sealed squeeze- film 


It is not uncommon to use seals in squeeze-film dampers to 
reduce end-leakage and thus to increase the effective damping. 
This application, however, makes computations difficult to 
carry out since any mathematical model for predicting the 
squeeze-film pressure is dependent on the amount of oil leakage 
allowed by the seals. First, let us consider that the seals 
allow only small leakage from the ends of the damper. It is 
therefore assumed that the pressure gradient in the axial 
direction is negligible compared to that in the circumferential 
direction. This is known as the long-bearing approximation. 
Hence eqn.(l) becomes 

9 _3 9 p 

g-g(h g-g) = 12 (e ' COS 0 + £ 01 ' S i 0 ) ( 3 ) 


The solution_ of_eqn ( 3 ) for pressure, satisfying boundary 
conditions p = p Q at 0 = 0, and 2ir is readily available in the 
literature [2] and is quoted here as 


P c (e) 


- 12 a' ( £ - — ) 

2 + £ 


(2 + £cos0)sin0 
2 

(1 + ECOS 0 ) 


+ 


6f'[ 2 

(1 + £ CO S 0 ) ' 


(1 + 


2 N 2 
e ) 


J + P, 


( 4 ) 


where p 0 is the nondimensional pressure at the maximum film 
thickness. It could be taken as the supply pressure p s but 
more reasonably can be obtained from the short-bearing pre ssure 
expression (eqn (2)) by setting the conditions p(0,z) = p Q at 
0 = 0 , and averaging over the land length. This gives 


- (i) 2 


(1 + e) 


J s' 
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arid ensures that p c (0) in eqn (4) is not a function of z. 

p s is used instead of j^/ 2 since with sealing there will be 
little attenuation of p s . Eqn (4) is true only for squeeze- 
film bearings in which the pressure gradient in the axial 
direction is relatively small, or for bearings which have a 
large length to diameter ratio. 

It appears that neither the short-film model nor the long- 
film model offers a suitable solution from which to calculate 
the pressure distribution in a centrally-grooved and sealed 
squeeze-film damper with a restricted but not necessarily small 
oil leakage. This is because the pressure gradients in both 
axial and circumferential directions should both be 
significant for such dampers. For this reason and for speed of 
computation using non-linear programs, an analytical pressure 
expression which includes the effect of axial oil flow as well 
as circumferential oil flow is essential in order to calculate 
the squeeze-film forces. 

To construct such an analytical model, let us assume that 
the pressure distribution in the z-direction in a centrally 
grooved and sealed squeeze-film damper is to be built up as 
illustrated in Fig.l. Now the boundary condition at z = - 1/2 
is readily available as the oil supply pressure in the groove, 
hut the boundary condition at the ends of the squeeze-film (at 
z = 1/2) is dependent on the oil-leakage flow rate. For a 
given damper design this boundary condition may be specified by 
introducing an end-leakage factor, say X, and writing the 
boundary condition as 


P = X p c (e) , 

in which X = 0 for full leakage and X = 1 for zero leakage. 

p c ( 0 ) is readily given by eqn (4) as the nondimensional 
pressure in the 0-direction for the relatively simple case of 
the long-bearing approximation. It should be noted that this 
allows pressure to be variable in the circumferential direction 
at the ends of the damper, whereas often the short-bearing 
approximation is used in conjunction with the specification of 
a constant pressure at the ends of the damper. 

Now using these boundary conditions as 

"p = "is at T - - 1/2 

and p = Xp c ( 0 ) at "z = 1/2 , 


we obtain 
"p(0/z) = 


the final form of the modified pressure expression 


>,L S 2 ( e 1 cos0 + ea'sinO) —2 1 N , , — x , — . l x 

6<s) ,i + ccose) 3 U ' + 1P = <9) <Z + 


+ p s (^ - z) 


(5) 
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Fig. 2(a) shows a typical set of circumferential pressure 
distributions at the mid-land position for two different values 
of the factor X, together with the theoretical predictions for 
the short-bearing and the long-bearing film assumptions. Fig. 
2 (b) shows a comparison of axial pressure distributions between 
the short-bearing and the modified film model for X = 0.09. 

The indication is that the modified film model may provide a 
useful approximation to the experimental case when the value of 
X is selected correctly. 

Calculation of squeeze-film forces with curtailment of 

negative pressure 

The squeeze- film forces and P 2 are obtained by 
integrating the squeeze-film pressure distribution in 
directions respectively along and normal to the line of centres 
of the journal and the bearing as follows: 

P l,2 = " ^ ( § )2rL / 0 " P (fl ' 7) (sine) dTd9 (6) 

where 0 i and 62 are the film boundaries in the circumferential 
direction . 

In this work, with a view to obtaining adequate 
predictions of vibration orbits, each numerical hydrodynamic 
pressure distribution was curtailed at its experimentally- 
observed negative pressure limit before integration to obtain 
the squeeze-film forces as in ref [1]. 

Experimental tests 

Tests were first carried out using oil of viscosity 
6.15 cp, an end-plate clearance of 0.0635 mm and a supply 
pressure of 5 lbf/in 2 gauge (345 kN/m 2 ) . Fig. 3(a) shows 
counter-clockwise orbits within the clearance circle of 0.216 
mm radius, and dynamic pressure recordings at a mid land 
position at the base of the damper, at various rotor speeds for 
Q c = 0.459, while Fig. 3(b) shows the corresponding numerical 
predictions, for which the value of the X factor used was 0.09. 
The negative pressure limit noted in each experimental pressure 
recording was used as a limiting condition for the corres- 
ponding numerical computation of dynamic pressure. Similar 
comparisons for a higher value of Q c of 1 . 055 are given in 
Figs. 4(a) and (b) again using a X factor of 0.09 for the 
numerical computations . These results show experimental 
vibration orbits and pressure waveforms in good agreement with 
their numerical counterparts. In all the experimental pressure 
waveforms the horizontal base line denotes atmospheric 
pressure . 

Effect of supply pressure on the damper performance 

Since the supply pressure is added linearly to the dynamic 
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variation in the X factor appears to be relatively small, in 
particular for tight sealing, that is low values of d/c. 


SOME FURTHER CONSIDERATIONS OF THE OUTLET BOUNDARY CONDITION 

As already pointed out the outlet boundary condition for 
the sealed squeeze-film is dependent on the oil leakage flow 
rate allowed by the seals, A means of specifying this boundary 
condition was given by the introduction of an end leakage 
factor X. Some further considerations relating to the 
introduction of this X factor will now be given. 

Assume a solution to the Reynolds equation (eqn 1) of the 

form 


p(e,z) = f(z) p c (6), (7) 

where ffz) and"p^-,(0) are functions for pressure in the axial 
and circumferential directions respectively, the latter being 
specified at the mid- land position. Substituting eqn (7) into 
eqn (1) and dividing by p c (0)h^, gives 

b" + g 2 ffz) - =0 (8) 

dz' 

where b = 12 (e ' cos0 + ea' sin0)/p c {0) TP* , and 


_ i_ 

3 0 (h 


0P C ( ® 

“a? — ) /p c (e) h 


The solution of eqn (8) for f (z) , with g 2 being treated as a 
parameter, is straight forward, and is of the form 

ffz) = Cl e gLT/E + C 2 ,- gLT / R - \ (9 

where C^ and Cj are constants to be determined from the 
boundary conditions. For the damper under consideration these 
boundary conditions are: 

7(0,^) = "p's at "z" = - 1/2, 

and a boundary condition at z = 1/2 which can be obtained by 
considering the axial flow balance across any unit 
circumferential area around the seal as 


where 


h 


3 


12n 




12 tU 


/ 


d is the end seal gap width, 

Z is the effective length of the gap. 
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pressure, the effect on the vibration orbits of a higher supply 
pressure should theoretically be akin to the effect of a lower 
negative pressure. In fact, an increase in the supply pressure 
would result in a reduction in the orbit size, and sufficient 
increase would suppress the negative pressure region completely 
and thus eliminate the load carrying capability of the squeeze- 
film [3] . To show the effect of supply pressure on the 
vibration orbits and cavitation pressure, a set of experimental 
recordings is given in Fig. 5(a) for Q c equal to 1.055. 

Fig. 5(b) shows the corresponding numerical predictions for 
X = 0.09. For these tests supply pressures of 5, 25 and 45 
psig (34.5, 172 and 310 kN/m 2 ) were used. From these figures, 
it can be seen that the reductions in the orbit size are 
substantial between 5 and 25 psig , whilst the reductions 
between 25 and 45 psig are relatively small . An accompanying 
feature is the development of negative pressure with increase 
in supply pressure, as can clearly be seen in the recordings of 
dynamic pressure. This may be due to the fact that for 5 psig, 
venting from atmosphere takes place, while for higher supply 
pressures this is prevented by continuous flushing. Such 
flushing would remove cavitation bubbles and allow the tension 
spikes to form as shown at the higher supply pressures. It is 
also interesting to note that no significant rise in the 
positive peak pressure level was observed or predicted with 
increase in supply pressure, probably due to the accompanying 
reduction in orbit size. 

Effect of the sealing gap on damper performance 

The next phase of the present work was to investigate the 
effect of end-plate sealing on damper performance and thus to 
give a graphical estimation for the value of the X factor for a 
given end-plate clearance. For convenience testing was first 
carried out for a supply pressure of 5 psig. (345 kN/m 2 ). Three 
different end -plate clearances of 0.114 mm, 0.139 mm and 
0.216 mm were used, being set by the insertion of spacing 
shims. The experimental recordings shown in Fig. 6(a) are for 
a Q c value of 0.549, a speed of 3500 rev/min and for an oil 
vicos ity of 6.15 cp. These compare with the display of Fig. 

3(a) (3500 rev/min) for a sealing gap of 0.0635 mm. Fig. 6(b) 
shows the corresponding numerical predictions. For each 
sealing gap the value of the most appropriate X factor is given 
in the figures. From these comparisons, it can be seen that 
the general size of the orbit increases with increase in the 
value of end-plate clearance. Good agreement can also be 
observed between the experimental recordings and the numerical 
predictions for both orbits and pressure waveforms . From a 
wide range of results, covering a variety of rotor speeds, 
supply pressures and unbalance ratios u/c,a plot of the X 
factor against the end-clearance ratio (defined as the ratio of 
end-plate clearance to the damper radial clearance) is shown in 
Fig. 7. It indicates that for a given end-clearance the 
appropriate value of the X factor is increased with increase in 
oil viscosity. However, given the viscosity range used, the 
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p and 3p/3z are the pressure and the pressure 
gradient at the end of the damper (z = L/2) 
respectively. 

This boundary condition has been used by Marmol and Vance [5] 
with some success . Rearranging and writing in nondimensional 
form, we obtain, 

p ( 9 , z ) = p = - yh 2 (4^) at z = l/2, 

3 z 

where y = (c/d) 3 Jt/L, and is introduced as a flow coefficient. 

Thus, in terms of the function f(z), these boundary 
conditions are : 

f(I) ='Fs/Pc(9) at ”z = - 1/2, and 

f(z) — - yh 3 df/dz at z = l/2. 

Inserting the above boundary conditions into eqn (9) and 
solving for Ci and C 2 , we obtain 

1L) e gL/2R - C.e gL/2R , (10 

g 2 


where 

(b ) e gL/2R _ L P s . + b -jd + yg Y 3 L /R) e 3 9L/2R 
9 P c ( e ) g 2 

(1 - yh 3 L/R) - (1 + yglT 2 L/R)e 2gIj ^ R 

Equation (9) for the f(z) function also requires the functional 
form of g 2 . Observing that as L/R -► <», df/dz 0 and f -»• 1 , 
then from eqn (8) g 2 becomes 

g 2 = - b = - 12 (e'cos9 + ea ' sin6 ) /p c ( 0 Jli 3 (11) 

where p c (6) is the nondimensional pressure from the long- 
bearing approximation of the Reynolds equation. Having 
obtained the required functions the pressure can then be found 
from eqn (7) as 

IpU/z) =~p c (e) (Ci eSWR + C 2 e-g L z/R + 1 ) (12) 

where g, Cj and C 2 are calculated from eqns (11) and (10), 
respectively. - Equation (12) is a general pressure expression 
for any finite- length squeeze-film damper with central 
circumferential oil supply and end-seals . A similar approach 
was used by Barrett et al [6] but for an open-ended journal 
bearing. It was shown to give reliable results. 





Pc( 0 ) 
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Some idea of the variation of the pressure distribution 
with sealing can be obtained from Figs. 8(a), (b) , and (c) . 

Fig. 8 (a) shows a typical set of circumferential pressure 
distributions at the mid-land position for two extreme values 
of the flow coefficient y (y = 0 and y -»• «>), together with the 

theoretical predictions based on the short-bearing and the 

long-bearing film assumptions. In this comparison the modified 
film model is also included, using a value of the X factor of 
0.09. Fig. 8(b) shows a comparison of axial pressure 
distributions for these models, whilst Fig. 8(c) shows a 
comparison of circumferential pressure distributions at the 
ends of the damper for the finite-length and the modified film 
models. It should be noted that the value of y of 27 

corresponds to the case of tight sealing for which the value of 

X of 0.09 was found by trial to give the best agreement in 
respect of orbit size and pressure wave form. From Fig. 8(a), 
it can be seen that even with the full sealing condition 
(y -»• «>) the corresponding peak pressure is considerably smaller 
than that obtained from the long-bearing approximation. For 
full leakage (y = o), the finite-length solution is nearly 
identical to the short-bearing solution, as is expected due to 
the small L/R ratio employed in the test bearing. As far as 
the specification of the end boundary condition is concerned, 
similarity (but not equality) can be observed in Fig. 8(c) 
between the pressure profiles obtained from the two different 
approaches of using the X factor of 0.09 and the y coefficient 
of 27. 

To give a graphic view of the approach, two sets of 
numerical orbits and pressures are shown in Figs. 9(a) and (b) , 
using respectively values of y of 27 and «> (the full sealing 
condition) . These results should be compared with the 
experimental recordings of Fig. 3(a) and the numerical 
predictions of Fig. 3(b). From these comparisons, it can be 
seen that, even with full-sealing (y = »)the predicted orbits 
are considerably larger than the experimental ones and the ones 
predicted by using the X factor of 0.09, in particular for the 
case of 5800 rev/min. A reversal in the position of the sharp 
"tail" between 3500 and 4000 rev/min was not observed when 
using the y coefficient. The numerical pressure waveforms of 
Figs 9(a) and (b) also do not show good agreement with 
experiment. These comparisons are typical of the other 
experimental cases considered. A general observation was that, 
using the approach of this appendix, the numerical predictions 
were in poor agreement with the experimental results, in 
particular for high values of Q c and at high rotor speeds, with 
tight sealing. 

The X factor discussed previously may be introduced into 
the specification of the outlet boundary condition by again 
considering the end pressure obtained from the flow balance as 
follows. This pressure has been shown to be, 

p( 0 , z ) = p = - yh 3 (4E) at z = 1/2 

3 z 
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From eqn (7) the pressure gradient is 


= pT(e)(4I>- 

3 z c 3 z 

Hence, 

P « p_ ( 6 ) ( - yh 3 41) (13) 

c 3z 

Comparing this with the equation, 

"p = Xp c <0)» (14) 

shows that A replaces 


(- Yh 3 41 ) 


There are two possible reasons why equation (13) does not give 
satisfactory results. These may be enumerated as follows: 

1) The leakage flow through the end seals may not be truly 
radial. However, even if a screw type of flow is assumed 
by increasing y say, we have seen (Fig. 9(b)) that this 
still does not give adequate agreement with the 
experimental results, even when y — 00 • 

2) Equation (13) assumes reverse flow from the seals 

when the pressure p is negative. However, this does 

not seem likely since there is no oil reservoir behind the 
seals, and indeed in the experiments, ingress of air was 
the more likely occurrence. 

Thus, while there is no obvious parametric link between A 
and the dimensions of the end seals, the A approach appears to 
be one to give adequate predictions of vibration orbit and 
dynamic squeeze-film pressure. Fig. 7 can then be used to 
establish an empirical link between parameters. 

CONCLUSIONS 


The work presented in this paper has dealt with the 
performance of the sealed squeeze-film damper when used in 
aero-engine applications. In such applications, the operation 
of the damper is governed by five nondimens ional groups. These 
groups are the nondimensional static force Q, the nondimen- 
sional dynamic force Q c , the nondimensional viscosity 3, the 
nondimensional supply pressure p s and the nondimensional 
cavitation pressure p cav . Good agreement has been observed 
between experimental observations and numerical predictions 
using an empirical leakage factor, over a wide range of 
operating parameters . This work has also demonstrated 
experimentally that sustained lift of the rotor in the squeeze 
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film annulus can be obtained without the use of any parallel 
support stiffness, and the squeeze-film was shown to provide 
this lift by the assumption of a cavitated film model. 

It has been shown that the use of an end-leakage factor X 
has resulted in successful predictions of vibration orbits and 
dynamic pressures in the sealed squeeze-film. It has been 
demonstrated experimentally that the damping can be increased 
effectively by decreasing the sealing gap and a means of 
determining X for a given gap has been given graphically. One 
indication was that this gap should be made less than half the 
damper radial clearance to achieve effective damping. 

There remains the problem of needing to know in advance 
the value of negative pressure to adopt for curtailing the 
squeeze-film pressure during computation of the squeeze-film 
forces and orbits in any given application. Experimental 
recordings of dynamic pressure show that this negative pressure 
appears to take one of three approximate values : 

(i) Approximately atmospheric, when the supply pressure 
is not great enough to prevent ingress of atmospheric 
air, which fills the cavitation zone; 

(ii) Approximately absolute zero, when cavitation bubbles 
exist due to the liberation of dissolved gases and 
fluid vapour. The supply pressure here has to be 
somewhat higher to prevent ingress of air, but is 
still not high enough to flush the bubbles away; 

(iii) Lower than absolute zero, when the oil temporarily 
supports tension, many cavitation bubbles in the 
previous cycles being flushed away by a fairly high 
supply pressure. 

For the bearing assembly considered in this work, when 
used under most engine conditions, it could reasonably be 
expected that much of the performance could be such as to 
produce cavitation pressures in the third region when a value 
of about -40 psig seems a fair average. This is based on the 
assumption that the external oil pressure at the squeeze-film 
groove is likely to be around 40 - 50 psig. It is suggested 
that, given the highly complicated mechanism of negative 
pressure development in the squeeze- film, reliance can be 
placed on the use of fairly high supply pressures to ensure a 
measure of consistency in the negative pressure at which film 
rupture is assumed to occur. 


418 



REFERENCES 

[1] Holmes, R. and Dogan, M. ' Investigation of squeeze film 
dampers in flexible support structures ' , NASA Conference 
Publication 2250 entitled 'Rotor dynamic instability 
problems in High Performance Turbomachinery - 1982 V. 

[2] Pinkus, 0. and Sternlicht, B, 'Theory of hydrodynamic 
lubrication' . McGraw-Hill Inc., New York, 1961. 

[3] Holmes, R. 'The non-linear performance of squeeze-film 
bearings'. Jnl Mech.Eng .Sci . , Vol . 14 , no.l, 1972. 

[4] Holmes, R. and Dede, M. 'Dynamic pressure determination in 
a squeeze-film damper'. Inst. Mech.Engrs . Second Int. 

Conf. on Vibrations in Rotating Machinery, Sept. 1980. 

[5] Marmol, R.A. and Vance, J.M. ' Squeeze- film damper 
charactersitics for gas turbine engines'. ASME Paper No. 
77-DET-18 , Oct. 1977. 

[6] Barrett, L,E,, Allaire, P.E. and Gunter, E.J. 'A finite 
length correction factor for short bearing theory ' . Trans . 
ASME Jnl Lub. Tech, Vol. 102, July 1980. 



Fig. 1 Damper configuration with circumferential 
oil supply groove and end plates 
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